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A Green-function formalism for the Kondo lattice model is presented, which is designed 
to be combined with the dynamical mean-field theory. With use of Wick's theorem only for 
conduction electrons, dynamical quantities are represented in terms of the t-matrix and its 
two-particle generalizations. By taking the high-frequency limit of the t-matrix with respect 
to a part of the fermion frequencies, one obtains dynamical correlation functions of localized 
moments. Explicit examples of calculational steps are provided with use of the continuous-time 
quantum Monte Carlo method. 
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1. Introduction 

In this paper, we present a formalism which is designed 
to be apphcable to the Kondo lattice and related mod- 
els with localized degrees of freedom. The Kondo lattice 
is the simplest model to investigate rich consequence of 
the collective Kondo effect such as heavy electrons, and 
competition with magnetically ordered states. In higher 
dimensions, the dynamical mean-field theory (DMFT) 
is the most effective and the simplest approach.^ If one 
tries to apply the DMFT to the Kondo lattice, however, 
a difficulty arises because the localized moments cannot 
be dealt with by the ordinary Green functions. This dif- 
ficulty results from the special situation of the strong 
correlation limit where the charge degree of freedom has 
been eliminated. Namely the fermionic Green function 
cannot be defined for the localized degrees of freedom. In 
this sense, the Anderson lattice is simpler since the Green 
function of localized electrons is well defined. However, 
the strong correlation limit of the Anderson lattice is nu- 
merically more difficult because the high-energy interme- 
diate states are essential in addition to singly occupied 
states for each site. Therefore, we pursue a formalism 
which allows one to perform highly accurate numerical 
calculation dealing only with low-energy localized states. 

In addressing to the localized limit, it is instructive to 
compare with the local Fermi liquid theory, which de- 
scribes low-temperature properties of the Kondo prob- 
lem. Nozieres has described the low-energy excitations 
in terms of the phase shift of conduction electrons near 
the Fermi level. ^ The phase shift contains all information 
of the impurity scattering caused by the exchange inter- 
action. By expanding the phase shift with respect to the 
quasi-particle distribution, one obtains the Wilson ratio 
as well as the transport properties at low temperatures. 
The same consequence can be derived by means of the 
perturbative approach of the Anderson model. Our ap- 
proach is a microscopic version of the phase-shift theory 
in the sense that we make maximum use of the viewpoint 
from conduction electrons. 



The impurity i-matrix corresponds to the product of 
the 4/-electron Green function and the hybridization 
squared in the Anderson model. Even the localized and 
strong-correlation limit keeps the t-matrix finite. Hence, 
with use of the impurity t-matrix, a unified treatment 
of single-particle properties is possible for the Ander- 
son and the Kondo- type models. In the same spirit, we 
will introduce a "generalized t-matrix" to describe two- 
particle correlations due to the local interactions. It will 
be demonstrated that the generalized ^-matrix gives the 
susceptibility of the localized moments both in the An- 
derson and Kondo-type models. 

The ideas of the impurity t-matrix and the general- 
ized t-matrix are also applicable to the periodic model. 
With the DMFT, the Kondo lattice model is mapped to 
the corresponding effective impurity model. The effec- 
tive bath is determined through the impurity t-matrix 
in a self-consistent manner. Furthermore, we will derive 
spatial correlations of the localized moments within the 
DMFT. The procedure which uses the Bethe-Salpeter 
equation with the local vertex of the effective impurity**' ^ 
will be extended to the Kondo lattice with use of the gen- 
eralized t-matrix. 

This paper is organized as follows. We present impurity 
and periodic models in the next section. In §3, the t- 
matrix is extended to the two-particle Green function 
with use of the equations of motions. The DMFT self- 
consistent equations are given for the localized models 
in §4. We further derive spatial correlations of localized 
moments within the DMFT. In §5, we present how to 
evaluate the generalized t-matrix in the continuous-time 
quantum Monte Carlo method (CT-QMC). An example 
of the generalized t-matrix is explicitly given in §6. We 
finally summarize in §7. 

2. Models 

The Anderson lattice model has been most widely used 
to discuss properties of 4/-electrons systems. We con- 
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sider a A''-fold degenerate band represented by 

Hc = J2^kci„(^ka, (1) 

ka 

where c^^ and Cka stic the creation and annihilation op- 
erators for the conduction electron with a-th compo- 
nent, respectively. The A^-fold degenerate Anderson lat- 
tice model is given by the following Hamiltonian: 

Fal = ffc + e/ E flfi- + ^ E(/^^- + 

ia ia 

U ^ ] fioifioificx'fia'i (2) 
i,{oia') 

where f}^ and fia creates and annihilates the localized 
electron at i-site, respectively. The sum over i is to be 
taken over all the 4/ sites, and {aa') denotes a pair of 
different components. 

By eliminating the charge degree of freedom, we de- 
rive effective models in the localized limit. We first ex- 
clude double and higher occupations by taking U = oo. 
Then taking the limit, V ^ oo and ef — oo with 
J = — V'^/cf fixed, the second order perturbation theory 
leads to the Coqblin-Schrieffer (CS) lattice model^ 

HcSL = i?c + ^ E fLf^c'4„'Cia, (3) 
iaa' 

which describes interaction between localized and con- 
duction electrons due to the virtual excitations to the 
4/° state. Since the intermediate 4/° state is isotropic, 
the interaction has SU(A'^) symmetry. If doubly occupied 
states are taken into account as the intermediate states, 
the interaction becomes anisotropic and more compli- 
cated. 

In the case of = 2, we can obtain a model with 
higher symmetry by including virtual excitations to dou- 
bly occupied state, since 4/^ state is isotropic as in 4/°. 
In the symmetric condition U = — 2e/, we take the local- 
ized limit, i.e., V — > oo and e/ — > — oo keeping the ratio 
J = —2V'^/ef. Then, we obtain the Kondo lattice model 

HKi. = H, + jJ2Si-al (4) 

i 

where Si is the localized spin at z-site and <t^ is defined 
by ctJ = X^^^/ cl^<T„„'Cia- The interaction term is invari- 
ant under the particle- hole transformation. The exchange 
interaction relates to the CS interaction as 

(7(7' <T 

where we have used the condition f^f^=\. The sec- 
ond term is the potential scattering independent of the 
spin component. The same potential for all sites are 
equivalent to a shift of the chemical potential, so that 
eq. (4) is rewritten as 

HkI. = E(^fe + ^)4<TCfea + JJ2 (6) 

where v = —J/2. The Kondo lattice model can there- 
fore be dealt with by replacing the chemical potential fj, 
with fj,' = 1^ — V in the CS lattice model. Hence there 
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is essentially no difference in the physical properties be- 
tween the Kondo and the CS interactions with N = 2. 
In this paper, we employ the CS interactions, which in- 
cludes N > 2. We note that only the t-matrix differs in 
the two models since they depend on the choice of the 
bare Green function. 

The above Hamiltonians of the lattice give the corre- 
sponding impurity models, if the index i is restricted to 
a single site. Hereafter, the terms 'Anderson model' and 
'CS model' refer to both the corresponding impurity and 
lattice models. We use 'impurity' or 'lattice' explicitly to 
specify either systems. 

3. Dynamics of Conduction Electrons and Gen- 
eralized f-matrix 

In this section, we present a formalism based on the 
Green function to handle the Anderson and CS models 
in a unified way. In the atomic limit, it is necessary for 
Green-function formalism to represent all quantities in 
terms of the conduction electrons. With use of equations 
of motions for conduction electrons, we derive dynamics 
of localized moments. 

In the finite-temperature formalism, it is convenient to 
work in the grand canonical ensemble. To this end, we 
replace the kinetic energy term He by 

He - M-^c = E ^kcl^Cka, (7) 

where ^fe = Cfc — /x is the kinetic energy with respect to 
the chemical potential. 

3. 1 Single-particle Green function 

We begin with the single-particle Green function prior 
to two-particle counterparts. With use of the Hciscn- 
berg operator Cka{T) = e^^Ckae~'^^ , the single-particle 
Green function is defined by'"'^ 

Glak'(.'ir,r') = -{TrCka{r)cl^,{T')), (8) 

where Tt is the time-ordering operator, and the bracket 
denotes the thermal average, (• • • ) = Tr{c''(^"-^) • • • } 
with e~^^ = Tii:e~^^ . In impurity systems, the site- 
diagonal element is given by 

G^a' (r, t') = N^' J2 Gkak'a' (t, t'), (9) 
fefe' 

where Nq is the number of sites. The Fourier transform 
is defined by 

(10) 

where e„ = (2n + l)7rT is the Matsubara frequency for 
fermions. We express the Green function of free conduc- 
tion electrons by gk{ien) = (ie„ - ^k)~^- 

To sec influences of the localized moment on the con- 
duction electrons, we derive equations of motion for 
^kak'a' ('''' ''"')• ^o this end, we first show equations of mo- 
tion for Cfe(r) and c^,(r). The Heisenberg equation leads 
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Fig. 1. Diagrammatical representation of the equation for 
G"^(ien)- Thick and thin lines denote the full and bare Green 
functions, G'^(ien) and g{ie„), respectively. 



to 



dcka (r) 
dr 

dr 



(11) 



where jfca is defined hy jka = [cfea, i?ioc] witli iJioc = H- 
He being the local interaction part of the Hamiltonian.^ 
Explicit forms of jka for each model introduced in §2 are 
given as follows: 



by 

G'fca(ie„) = .gfca(ie„) + gfcQ(ie„)tfea(ien).gfcQ(ie«), (15) 

tka{r,T') = -{Trjkc.{r)jL{T'))+S{r-r'){{jka,cij). 

(16) 

In the Anderson lattice model, the i-matrix corre- 
sponds to tkai^^n) — ^^^'{^.(ien) as in the impurity 
model. In the CS lattice model, the second term gives 
(J/Nq) J2k'{fk'afk'a)- By means of the equations of mo- 
tion, an expression of the internal energy can be derived 
from the single-particle Green function (see Appendix 
B). 

Finally, we conclude by noting the relation between the 
t-matrix and the occupation number of the local states. 
Since the t-matrix in the Anderson model is identical to 
the /-electron Green function with the factor V'^, the 
/-electron number for each component, n^, is given by 



Jfco 



N, 



-1/2 



Vfa, (impurity Anderson model) T ^ t„(ie„)e'^"' = F^^^^ (Anderson model) (17) 



Vfka, (Anderson lattice model, eq. (2)) 



Nr 



k' qa' 



-1/2 



J ^ fljaCa', (impurity CS model) 



k' +qa^h~qa' ■ 

(GS lattice model, eq. (3)) 



(12) 



By using jka, all models can be treated concurrently. 
Since the Kondo exchange interaction can be expressed in 
terms of the GS interaction as described in §2, we employ 
the GS model for definiteness. For convenient description 
of the impurity models, we introduce the impurity-site el- 
ement ja — X^fc Since the interaction is local, 

jka is independent of k in the impurity models and j^ 
can be simply written as ja = V Ngjka ■ 

Generating an equation of motion for G'f,^i^,^,{T,T') 
and taking the Fourier transform with respect to r — r', 
we obtain an expression of G^^^,^,(ie„). For impurity 
models, the site-diagonal component, G^(ie„), is written 
in terms of the impurity i- matrix ta(ien) as 

G^(ien) = gQ(ie«) + ffQ(ien)ia(ien)ffa(ien), (13) 

t„(r,r') = -(r.Ja(r)jl(r')) +<5(r-r')(0„,4}). 

(14) 

Its derivation is given in Appendix A. Figure 1 shows 
the corresponding diagrammatical representation. The 
impurity i-matrix includes all consequences of impurity 
scattering. The argument ofta{e+iO) — |iQ(e)|e"^'^'^) gives 
the phase shift of the conduction electrons (see for exam- 
ple ref. 10). In the impurity Anderson model, the second 
term vanishes and the t-matrix is given by the /-electron 
Green function, taii^n) — y^G^(ie„). In the impurity CS 
model, on the other hand, the second term gives the con- 
tribution in the first Born approximation J(/j^/a), which 
is independent of energy in proportion to an occupation 
number of the state a. 

For the periodic model, the Green function becomes 
diagonal in the momentum space, G^„(ie„), and is given 



where 5 is a positive infinitesimal. In the CS model, on 
the other hand, the Born term refers to the occupation 
number. Since the first term in eqs. (14) and (16) van- 
ish in the high-frequency limit, the Born term can be 
extracted by 



lim ta(ie„) = Jn^. (GS model) 



(18) 



These ways to derive the occupation number from the 
t-matrix will be extended to the two-particle Green func- 
tion in evaluating the local susceptibility. 

3. 2 Two-particle Green function 

By writing down the equations of motion for the single- 
particle Green function, wc have elicited information on 
the local moment. Here we extend this procedure to 
the two-particle correlation functions. We define gener- 
alized susceptibility for conduction electrons with four 
time variables as follows: 

Xl234(n,T2,T3,r4) = {Trc\{Tl)c2{T2)cl{T3)c4{T4)) 

- (T,cI(ti)c2(t2)>(T,4(t3)c4(t4)), (19) 

where the indices symbolically represent the wavenum- 
ber and the spin component. The Fourier transform is 
defined by 



xW(ien,ie 



1 

= In dri • • • / dT4Xi234(Ti, 7-2, ra, T4) 
P Jo Jo 

X eie.(r2-ri)gie„,(r4-r3)gii.„.(r2-T3)^ (2O) 

where — 2m7rT is the Matsubara frequency for 
bosons. The susceptibility is obtained by letting ti ~ T2 
and T3 = T4 in eq. (19), so that its Fourier transform 
Xi234(i^ni) is evaluated by 



X1234 



(ii/,„) = T^x^234(ienIien';i^^m)• 



(21) 



By analogy with the single-particle Green function, 
we consider extracting information of the localized elec- 
trons from x'^- For this purpose, we define a generalized 
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lEy, lEj,' lEj, ie,j' 



ie„+iVm iE„+iv,„ 



lE,i+lV,,, 1E„+1V,,, 



Fig. 2. Diagrammatical representation of the equation for 
Xi234(i«".i«n';i'^m) in eq. (22). 



t-matrix 71234 (icn, icn'; i'^m) for the two-particle Green 
function as follows: 

Xl234(i^ni iCri'! i'-'m) 

= -<5,m'[^i4'5235i(ien).92(ien + ii^m) 

+ 5i4gi(ien) • 53(ien + i^'m)i32(ie« + ii^m),g2(ie« + ivm) 

+ .g4(ie„)t4l(ien)gi(ien) • '523.92(ien + Wm)] 

+ 5i(ie„)g2(ien + ii^m)g3(ien' + ii^m)54(ieri') 

X Ti234{ien,ie„']iv„i), (22) 

which is diagrammatically expressed in Fig. 2. The gen- 
eralized ^-matrix includes all effective interactions be- 
tween two conduction electrons via the localized elec- 
tron. From an equation of motion for x'^, we can obtain 
an explicit expression of the generalized i-matrix (see 
Appendix A for detail). In the imaginary-time domain, 
'?'i234(ien,ie«';ii'm) is given by 

'?'l234('''l,T2,r3,T4) 

(7Vjj(ri)j2(T2)i3(r3)j4(T4)) -t2l{T2,Ti)U3{T4,T3) 

+ Sin T2)S{t3 n){TAjl{T,), C2(ri)}{jt(r3), C4(r3)}) 

- 6{n - t,)S{t2 - T3){Tr{jl{n), C4(Ti)}{jt(r3), C2(r3)}) 

+ <5(ri - T2){Tr{jl{n),C2{T,)}jl{TMn)) 

+ 6{t3 - n){Trjl{n)j2{r2){jl{T3), C4(r3)}) 

- 6{n - T,){Tr{jl{n),c,{n)}jl{Ts)Mr2)) 

- S{t2 - Ts){Trjl{n)u{T^){jl{Ts), C2{ts)}) . (23) 

In this derivation, we have used the relations {jj,C2} = 
{c\,j2} and [ci,{jl,C3}] = [c{,{jl,C3}] = 0. The latter 
relation holds in the case where i?ioc is bilinear in c' and 
c. 

In the Anderson model, terms except for the first two 
vanish in eq. (23) because of the relation {jj,C2} = 0, 
which follows from eq. (12). Hence T corresponds to the 
two-particle Green function of the localized electron as 
'?i234(ien,ie«';ii'm) = 

X 1234(1^1 1 ifn'! ii^m)j where X1234 
is defined in a manner similar to Xi234- Therefore, in the 
Anderson model, the susceptibility X'^i^'^m) is evaluated 

by 

Xl234(i^m) = -^X!'^1234(ien,ie„';il'm)- 



(Anderson model) 



(24) 



In the CS model, on the other hand, all terms having 
the delta-functions contribute to the generalized t-matrix 
in eq. (23). Figure 3 shows the corresponding diagrams 



+ 



iEj,-i-iv,„ ie„'-i-iv,„ 




ie,,'-i-iv„, iE„+iv, 




+ 




1E„ IE, 



+ 



ie,,+iv,,, iE,,'+iv, 




Fig. 3. Diagrammatical representations of terms in the general- 
ized f-matrix, eq. (23). Wavy lines indicate the delta-function for 
a time variable. 



in the frequency domain. Correlation functions vanish 
in the high-frequency limit, while the Fourier transform 
of 6{t) is constant in the frequency domain. Thus it is 
possible to extract each term in eq. (23) by taking a limit 
for the corresponding variables. For example, the third 
term, which has S{ti — T2)S{t3 — t^), can be picked out 
by a limit with respect to e„ and e„' such as 



lim lim Ti234(ien, ie,. 



■*4-oo e / — + — 00 



/9 



dT2 







dT3 



"'0 
{flf2){flU 



(r./t(r2)/2(r2)/](r3)/4(r3)) 



^ii/„(r2-T3) 



Tj\(234{^i^m), (CS model) 



(25) 



which is nothing but the susceptibility of the localized 
electron. The second term in the middle row comes from 
the second term in eq. (23), whose high-frequency limit 
is given in eq. (18). We note that, in eq. (25), the signs 
of e„ and e„i should be different from each other. If we 
take another limit keeping the difference = fn' — ^n, 
the fourth term in eq. (23) remains, so that we obtain 

hm lim Ti234( 

€-n — >- + 00 e / — » + 00 



X{234(i^'m) - Xf432(ii'™') " Ulfi) ifsh) ■ 



(CS model) (26) 



In conclusion, in order to obtain the susceptibility of 
the localized electrons, we may evaluate the generalized 
t-matrix T. Eliminating the fermion frequencies by cer- 
tain ways, T yields the two-particle correlations of the 
localized electrons. The ways to obtain the susceptibil- 
ities are different in the Anderson and the CS models. 
This difference originates in the relations between the 
occupation number and the ^-matrix, eqs. (17) and (18). 

We set up an integral equation for 71234 (icn, icn'; ii^m)- 
It has been shown in eq. (24) for the Anderson model 
that T1234 corresponds to the generalized susceptibility 
of / electrons, Xi234(i^ni i^"'! i'^m), where the irreducible 
vertex part is defined by the Bethe-Salpeter equation. In 
a similar manner, we introduce an irreducible vertex part 
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ie„+iv„, ie„'+iv,„ 




r 



ie„"+iv„, iE„'+iv,„ 



Fig. 4. An integral equation for the generalized f:-matrix 

n 1 



of Ti234 by 

; i^'m)/l'2'3'4' (ien, iCn"; if^m) 

n" 1'2'3'4' 

X 7;/3/34(ie„",ie„';ii/m), (27) 
where 7^234 (ie„; ii^m) is defined by 

^I°234(ien;i^m) = "^41 (iCn )i23 (ie« +il^m). (28) 

Tlie above equation is diagrammatically represented in 
Fig. 4. In the Anderson model, T° and / correspond to 
y4^/o g^j^jj r-'^/V^, respectively, with T-f being the irre- 
ducible vertex part. 

With these formulae, we now describe the following 
expression of the susceptibility in the lattice: 

= (TV4,a('ri)cfc+q,a(T2)cI'+q,a'(^3)Cfc',a'(T4)) 

- Sqfi{Trci jTi)ckAT2)){Trcl,^^,{T3)Ck'^a'{T4)), (29) 

The Fourier transform, defined in eq. (20), leads to 
). To discuss physical responses to 
external fields, we introduce two response functions con- 
structed by linear combinations of spin components. One 
is the charge response function defined from the 

charge density operator J2a '^a ' ^^'^ another is the mag- 
netic response function ^'^'^'^s defined from the magnetic 
moment, man"^ with J2a ™" = 0- 

Xkk'qi^>^n:ien'-^J^m) = J2 ^kk' q,ao,' i^'^n, i<^n' -^l^m) , 

act' 

(30) 

c.mag/. ■ • \ \ c /■ • • \ 

Xfcfc'q (ien,ien';i'^m) = 2^m„m„'Xfefe'q,aa'(ien,ie„';i'^m)- 

act' 

(31) 

We assume the SU(iV) symmetry. Then the above quan- 
tities may be rewritten as 

Xkl^l = NXk^7 + - (32) 

x^;:^:t=^(^^i^^7-^Z)^ (33) 

where ^(^'^^'^^ and denote respectively the diago- 

nal and off-diagonal element of Xaa' spin indices, and 
Cn — N~^J2a''^a is the Curie constant. The physical 
response is evaluated from the generalized susceptibility 

by 



where 7 = chg, mag represents the charge and magnetic 
channels. 

We introduce the generalized i-matrix 
7^^./g(ie„, ien'; i^'m) for each channel 7 to describe 
the response functions of the localized moments. From 
eq- (22), Xkk'q gives TJ^,, as follows: 

Xfefc'q(ien>ien';ii^m) 

= -NCjiSnn'Skk'[9(k)g{k + q) 

+ g{k) ■ g{k + qft{k + q) + g{kft{k) ■ g{k + q)] 

+ 9{k)g{k + q)g{k' + g)g(fc')Tj^,q(ie„, ie„'; ii^„), (35) 

where four-dimcnsional vectors are introduced by fc = 
(A;,ie„), k' = (fc',ie„') and q = {q,ium). The factor Cj^ 
corresponds to Cat for the magnetic channel and to unity 
for the charge channel. In the Anderson model, eq. (24) 
yields the susceptibility of the localized electrons, so that 
we obtain 



f''^{w,n) = -^^Tq^(ie„,ie„';ii/™), 



y4 



(Anderson model) (36) 



where is defined by N^^ J2kk' '^kk'q- Here, we have 
assumed that the matrix elements of the magnetic mo- 
ment of the localized electron are the same as those 
of conduction electrons. In the CS model, on the other 
hand, from eq. (25), gives Xq'''(iJ^m) in the limit 

f 1 
Xq (ii'm ) = 77772 lim lim V (^^" ' i^"' ' i*^™ ) • 

(CS model) (37) 

By definition, the magnetic channel is zero in the CS 
model, so that we represent Xq"^'^^ simply by Xq here- 
after. 

4. Dynamical Mean-Field Theory for Localized 
Models 

In infinite dimensions, the self-energy part due to local 
interactions is also local. Accordingly, the periodic sys- 
tem with local interactions can be reduced to a single- 
impurity model in an effective medium. The DMFT de- 
termines the effective medium from the local self-energy 
in a self-consistent manner.^' Spatial dependences of 
two-particle correlation functions in the periodic system 
can be evaluated from the local vertices in the effective 
impurity system. 

In this section, we apply the DMFT formalism to the 
CS lattice, where the ordinary Green function cannot be 
defined for localized electrons. For this application, the t- 
matrix and its two-particle generalization, introduced in 
the previous section, play a key role. The self-consistent 
equation for the conduction electrons will be constructed 
with the i-matrix. Later in this section, we further derive 
spatial dependences of two-particle correlation functions 
with use of the generalized i-matrix. 



(34) 



kk' nn' 
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4-1 Self- consistent equations 

It is convenient to begin with the Anderson lattice 
model, since its self-consistent equation has been well es- 
tablished. In the DMFT, an effect of U at the surround- 
ing sites is incorporated in an unperturbed Green func- 
tion. It is referred to as a cavity Green function Q-^'^\z). 
Letting Q^{z) be the full Green function in the effective 
impurity system, the self-energy part T^-^{z) is given by 

Ef{z) = gf°izr^^gJ'iz)-\ (38) 

In the infinite-dimensional limit, the self-energy part be- 
comes local, so that ^-^{z) is equivalent to that of the 
periodic system. The site-diagonal Green function G^{z) 
of the periodic system is given in terms of S'^(z) by 



Gf{z) = j dep(e 



-ef-^f{z)~ 



z — e- 



(39) 

where p{e) = Nq^ i5(e — efc) is the density of states 
of conduction electrons. The self-consistency condition 
g^{z) = [z) renews cavity Green function as 

==G-^(z)-i + I]/(z). (40) 

This equation corresponds to a procedure which removes 
^^{z) at one site keeping ^^{z) at the surrounding sites. 
A self-consistent solution of the above equations yields 
the Green function of the lattice in the DMFT. 

In order to employ a perturbation theory from the 
atomic limit as an impurity solver, it is convenient to 
rewrite the above equations in terms of the Green func- 
tion of conduction electrons. This transform enables us to 
apply the DMFT formalism to the localized model, such 
as the Kondo lattice and the CS lattice. To this end, we 
introduce a cavity field g^^{z) for conduction electrons 
defined by ^■''"(z) = [z-ej-V'^g'=^{z)]-'^ .Thcscli-cncigy 
part is local in the infinite-dimensional limit, so that the 
site-diagonal element G'^{z) of the conduction-electron 
Green function is given by 



G 



\z) = J dep{ 



e) [^-e + /.-E^(^)]- 



The self-energy part ^^{z) is defined by 

t{z) 



E%z) 



(41) 



(42) 



z-ef-^f{z) l + g^O{z)t{zy 

where the impurity t-matrix t(z) is defined by t{z) = 
V^g^ {z), and is related to the site-diagonal Green func- 
tion g'^{z) in the effective impurity system by 

g%z) = 5^0(z) + 5^0(z)t(z)g="(z). (43) 

It is obvious from eq. (42) that S'^ (z) characterizes an im- 
purity scattering without repetition, while the t-matrix 
incorporates all scattering processes including repetitions 
of S'^(z). The self-consistency condition for the localized 
electrons, g-^ [z) = G^ {z), ensures a similar relation for 
the conduction electrons, 5'^(z) = &{z). Hence the cav- 
ity field ^'^'^(2;) is given from G'^{z) as 

a"°(^)-' = G'^(2)-' + S^(2). 



(44) 



Equations (41)-(44), written in terms of conduction elec- 
trons and the t-matrix, are identical with eqs. (38)-(40), 



written in terms of the localized electrons. The DMFT 
equation for the CS lattice model eventually becomes the 
same as that of the Hubbard model. The sources of the 
self-energy 1^ are, however, different: in the CS lattice 
model, E'^ signifies the impurity scattering and accord- 
ingly, is obtained from the impurity t-matrix in eq. (42), 
while in the Hubbard model the self-energy arises from 
the Coulomb interaction between conduction electrons. 
From eq. (39), the f-matrix of the lattice, t{z), defined 
by i{z) = V^G^iz) is related to G^iz) as 

i{z) = i:%z) + i:%zfG%z). (45) 

The impurity f-matrix t{z) equals to t{z) from the self- 
consistency condition. In the localized model, it is possi- 
ble to study formation of heavy quasi-particles through 
the t-matrix of the lattice. 

The DMFT imposes self-consistency on the site- 
diagonal element of the Green function, so that only fc- 
averagcd quantities appear in the equations. However, it 
is possible to see the fe-dependence through the integrand 
of the average. The Green function and the t-matrix for 
a given k are evaluated from T,''{z) by 



Gl{z) = [z-ek + ^i-E%z)]-' 



9k{z) 



l-5fe(z)E-(z) 



(46) 



E^(z) 



tk{z) = E=(2) + E^(2)^G^,(2) = , ^ (47) 

1 - 9k[Z)l^^[Z) 

where the fc-dependence enters only through the energy 
Ck of the conduction electrons. 

4-2 Spatial dependence of two-particle correlations 

Let us consider the two-particle Green function within 
the DMFT. Applying eq. (22) to the effective impurity 
system, the local component of the generalized t-matrix 
7^2c is defined by 

+ a'^°(ie„) • a"°(ie„ + lUmfiiien + Wm) 

+ Q''%i€r^)Hiien)■g^%ien+Wm)] 
+ a^°(ie„)a^°(ie„ + iiym)G'=° {iCn' + Wm)G'\ien') 

xT[2ci^en,ien';iiym)- (48) 

In a manner similar to eqs. (36) and (37), 7^^^ gives the 
local susceptibilities for each model. Noting eq. (43), the 
above equation may be rewritten as 

-|-7r(ie„; ii'm)[7^|^^(ie„, ie„/; iz^^) 

-^n«'^o;''(ie„; ii/m)]7r(ie„'; ii^m), (49) 

where we have introduced 

7r(ie„;ii/m) = -t/''°(ie„)tj''"(ie„ +iv„i), 
nioc(ie„; ivrn) = -G"'(ie„)G"=(ie„ -I- ii/^), (50) 

and, x^oc'^ is defined by x^oc^ = iVC^Hioc- 
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We direct our attention to the generalized t-niatrix in 
the effective impurity as well as the lattice systems. We 
set up an integral equation for T^^. Equation (27) holds 
for each channel independently, so that we obtain 



+ T ^ 7^°;'' (ie„; ii/,„)/''(ie„, ie„" ; \Vm)TZS^'^n" , ifn' ; i^^m) 

n" 

(51) 

where T^^'^ {ien \ it'™) is defined by 

^oc'^(i£n; i'^™) = -A^C^i(ie„)i(ie„ + Wm)- (52) 
In the Anderson model, relates to the irreducible ver- 
tex part r-'^^^ by r = Vf^^ /V^, and T" to the general- 
ized susceptibility without the vertex correction, X''^"'^; 
by T^''^ = V'^x^^'^ ■ infinite dimensions, only diagrams 
for intersite processes are those connected by two prop- 
agators. It follows that I'^ is a local quantity.^' There- 
fore connecting at each site by a pair of the lattice 
t-matrix leads to the generalized t-matrix of the pe- 
riodic systems. Taking the Fourier transform over the 
spatial indices, we obtain an equation for as 

77(ie„,ie„/;ii^„) = (5„„/T°'''(ie„; ii^™) 



7q"''''(ie„; ii',„)/'''(ie„, itn" ; WmjT^ 



where T^'"^ is defined by 
T°''^(ie„;ii'm) = 



(53) 



No 



^ifc(ie«)ife-i-q(ien + mn)- 

k 

(54) 

We represent eqs. (51) and (53) in a matrix form for the 
fermion frequencies as follows: 

^7 '7-0,7 



- T 

^loc ^ -'loc 



_ 7-0,7 777-7 



7-7 _ 7-0,7 7-0,7777-7 



(55) 



Eliminating I'^ from the above two equations, we obtain 
the following equation for T^: 

IVr' = KJ"' - KcT' + [T,"^'']-', (56) 
which holds for each value of the energy transfer inde- 
pendently. This equation is an extension of those in the 
Hubbard model'* and the Anderson lattice model^ to the 
CS lattice model. The corresponding equation for the An- 
derson lattice model has also been derived within the ex- 
tended non-crossing approximation (XNCA), where the 
irreducible vertex is explicitly given in the framework of 
the NCA.^'^ We again note that, in the CS lattice model, 
a limit of 7^™^*; gives the magnetic susceptibility of the 
localized moments in eq. (37). Although eq. (56) is an 
infinite-sized matrix equation, we actually do not need in 
numerical calculations to compute T in the whole range 
of frequencies. The matrix equation as well as the lim- 
iting operation in eq. (37) can be evaluated with small 
matrices efficiently (see Appendix C). 

We return our attention to the susceptibility of con- 
duction electrons. The susceptibility makes it possi- 
ble to discuss instabilities against the charge-density 



+ 



1E,,+IV„, 1£,,'+1V,, 



ie„-+iv„, iE,,'+iv,„ 



Fig. 5. An integral equation for the generalized susceptibility of 
conduction electrons, x'^{\tn, '^^n'\ ''^^in)- 



wave (CDW) and spin-density wave (SDW) states, cor- 
responding to the charge and magnetic channel, re- 
spectively. For the generalized t-matrix, we have intro- 
duced the irreducible vertex /''' to derive spatially de- 
pendent function. To obtain an integral equation for the 
conduction-electron susceptibility, we substitute eq. (51) 
into eq. (49). After some manipulations, we obtain the 
following equation: 



T^6nn' X\oc (ien ; ii^m )r''''' (ie„ , ie 



IXioc(iei 



(57) 

For more details and the explicit expression of V^'"^ , see 
Appendix D. Figure 5 represents the above equation dia- 
grammatically. Since F'^'^ consists only of local quantities 
in infinite dimensions, F'^'''' is also local. Hence, with use 
of the same vertex, Xq'^ can be constructed as follows: 

Xq'*'(ie„,ie„';ii',„) = (J„„'Xq"^^(ie«; ii'm) 



)F'''''(ie„, ie„" ; ii/„i)Xq''(ien", 



(58) 



iVC^Hq with 



where Xq"'^ is defined by Xq"''' 
Hq(ie„; ii/,„) = --^ ^ Gfc(ie„)G'fc+q(ie,i + ii^m)- (59) 



Eliminating F'^^'^ from eqs. (57) and (58), we obtain the 
following matrix equation for Xq^ 

ix^r' - ixioir' - [x?o"c"]"' + ixf^T'- (60) 

This equation is identical to that in the Hubbard model."* 
This follows from the fact that all the spatial depen- 
dence can be represented by G^, and therefore by Hg, 
in both the Anderson lattice and the CS lattice. Equa- 
tions (57) and (58) represent effective onsite interactions 
between conduction electrons resulting from the interac- 
tion with the local moment. The total susceptibility in- 
cluding {cli^afai fa') Can also be derived in an analogous 
way.^ 

We should note that the DMFT guarantees self- 
consistency only in the single-particle Green function. 
Consequently the two-particle correlation function does 
not satisfy self-consistency between the effective impu- 
rity and the lattice models. Namely, eq. (60) gives 



Xlo 



oc 7^ ^ Xq) 



(61) 



as in the ordinary random phase approximation (RPA) 
theory. Hence, the present formalism corresponds to the 
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RPA level in momentum dependence, but proper local 
dynamics is incorporated through the irreducible vertex. 

5. Evaluation of the Generalized t-matrix in the 
CT-QMC 

Recently, a new impurity solver called CT-QMC has 
been developed. ^"^"^^ The single- and two-particle corre- 
lations in the impurity CS model can be obtained ac- 
curately within statistical errors. In this section, we 
describe how to evaluate the generalized t-matrix of the 
impurity CS model in the CT-QMC. It will be shown 
that the high-frequency limit for Green functions can be 
taken strictly in this formula. 

To establish notations, we begin by summarizing the 
formula of the single-particle Green function. In the 
Monte Carlo simulations, a snap shot is expressed by an 
integer k and two sets of variables, (n, • • • , r^) for time 
and (ai, • • • , afc) for spin, which we symbolically repre- 
sent by /C hereafter. For a given configuration of /C, the 
Green function is given by 

Ga(T,r';/C) = ga{T - r') 

- ^9c,{t - Tj){M„)jig^{Ti - t'), (62) 

where the matrix Ma is defined by Ma = and 
{Da)ij = g{T^' - rj) with {r/'} and {rj} being cer- 
tain sets of t's in JC. The average over a Monte Carlo 
ensemble gives the physical Green function: G(r, r') = 

(G(r,T';/C))MC. 

We consider the generalized susceptibility for conduc- 
tion electrons defined by 

Xaa'(Tl.T2,T3,T4) = (TVcl )Ca (r2)c^, (r3)c„/ (r4)) 

- {Trci{r,)Ca{r2)){Trcl{T3)Ca'{n)). (63) 

For each configuration of fC, the Wick's theorem is ap- 
plicable to the four-operator average, so that we obtain 

(^r 4 (n ) Ca (r2 ) 4/ (rs ) Cc,/ (r4 ) ; /C) 

= Ga{T2,Ti;IC)Ga'iri,Ts;IC) 

'Ga{T4,n;lC)Ga{T2,T3;lC). (64) 

This equation can actually be demonstrated with use of 

the fast-update formula for the four-operator addition.^** 
Performing the Fourier transform defined by eq. (20) 
and taking the Monte Carlo average, we obtain a for- 
mula for x^i^, (ie„, ie„/; it'm). The generalized f-matrix 
Xia' (ien, icn' ; i^'m), defined in eq. (22), is then given by 

Taa' (i^n ; i^n' ^ i^m) — {^a (^n : ^ ^m)^^a' i,^n' ~1~ ; ^n' ) 
- SmOta{^en)ta'{i£n'), (65) 

where 

Ua{ei,e2) = T^(M„),,e-^-^--i-'. (66) 

ij 

The impurity f-matrix is given in terms of Ua by 
ta(ien) = -(wa(e„,e„))MC- 



The generalized t-matrix. given in the frequency do- 
main, can also be expressed in the imaginary-time do- 
main. Prom eq. (65), we obtain 

Taa'{Tl,T2,T3,Ti) 



[{Ma)ji{Ma') {Ma)jl{Ma 




X 5{ti - Ti)5{T2 - Tj)6{T3 - Ti)6{t4 - Tm) 



ta{T2 - Ti)ta'{T4 - T3). 



MC 

(67) 



Although it is in principle possible to measure T in the 
imaginary-time domain, it is not practical due to its com- 
plex structure. Namely, two discontinuities for each vari- 
able between and /? make a precise treatment difficult. 
Instead, eq. (67) can be used for evaluating each term in 
eq. (23) selectively. For example, by letting ti = T2 and 
T3 = T4 in eq. (67), the third term of eq. (23) is singled 
out to give 

{Trfl{T)fa{T)flfa') 



J2 



J2 [iMc)ji{Ma>)ml 



- Sac'{Ma)jl{Ma)mi]S+{T,Ti - Tl) j , (68) 

/ MC 

where the /3-periodicity is imposed on S+{t, t'): 

0+[r,T)-y (t'<0) 

The summations in eq. (68) arc taken for all pairs satis- 
fying Ti = Tj and Tl = Tm- On the other hand, by letting 
Tl = T4 and T2 = T3, eq. (67) yields the transverse sus- 
ceptibility as 

{Trfl{r)fa'{r)flfa) 



T 



ral 



- ^aa'{Ma)jl{Ma)mi\^^{T, Ti - Tl) 



(70) 



MC 



As presented in rcf. 16, the susceptibility can also be eval- 
uated from the configurations, {t^} and {a^}, without the 
matrix Ma- While eq. (68) gives the same results as the 
configuration measurement within error bars, eq. (68) 
tends to have more errors. As for the transverse suscep- 
tibility, eq. (70) is the only way of evaluation. 

In order to obtain the spin correlations in the DMFT, 
we need the high-frequency limits of Taa'i}^ni ''^^n''i '^^m) 
as shown in §4. The high-frequency limit can be taken 
strictly in the same way as eqs. (68) and (70). In the 
limit e„ oo, only terms with 5(ti — T2) remain. Hence, 
this limit corresponds to restricting the summation in 
eq. (67) to Ti = Tj. After the Fourier transform, we 
obtain lim^^^oo '^a' (i^".> i^w; i^^m)- In a similar manner, 
the limit e„' — > 00 is calculated by restricting the sum- 
mation to Tl = Tm- 
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Fig. 6. 7^™'''^(ie„, ie„/ ; 0) for J = 0.3 and T = 0.02. 



6. Structure of the Generalized t-matrix 

In this section, we show a typical example leading to 
the q-dependent static susceptibilities. We employ the 
Kondo lattice model in the form of eq. (6). Namely, 
the bare Green function includes the potential scatter- 
ing V as gk{ien) = (i^n — Cfc — v). We note that, in this 
definition, the f-matrix does not have the particle-hole 
symmetry even at the half filling. We adopt the nearest- 
neighbor tight-binding band of the infinite-dimensional 
hyper-cubic lattice. "'^^ The density of states is given by 
p{aj) = D~^-\/2/7rexp(— 2a;^/D^), and we take D = 1 as 
the unit of energy. 

In the effective impurity problem, wc evaluate the 
static component, ivm = 0, of the generalized t-matrix of 
each channel 7, 7^^{i€n,i€n'',0). Figure 6 shows the mag- 
netic channel T^^^^'^iie,,, ie„'; 0) for J = 0.3 and T = 0.02 
at half filling. We plot only the range n' > because of 
the relation 7^^^(-ie„, -ie„'; 0) = T{^^{ien,i€nr,0)* . The 
frequencies arc taken up to n = 40 to see the whole 
structure of the function. As shown later, much fewer el- 
ements are actually sufficient for the evaluation of Xq- 
We notice, in Fig. 6, distinct structures around the diag- 
onal, n = n' , and the zero frequencies, n = or n' = 0. 
They are related to the ^-functions in eq. (23), which do 
not disappear even in the high-frequency limit. The fre- 
quencies can be exactly taken to infinity in the CT-QMC. 
Figure 8 shows 7ioc(ien) — ioo; 0) for the same parameters 
as above. While the function has some structure around 
n = 0, at high frequencies we confirm a convergence to 
Txioc, which is provided in eq. (25). On the other hand, 
the charge channel 7^^|!^(ie„, ie„'; 0) converges to zero at 
high frequencies, except for around the diagonal, n = n'. 
This high-frequency behavior follows because there is no 
charge degree of freedom in the localized electrons in the 
Kondo lattice model. 

In order to compute q-dependences of the generalized 
f-matrix with eq. (56), we need the lowest order quanti- 
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n 

Fig. 8. liocCien, -ioo; 0) for J = 0.3 and T = 0.02. 

ties without the local vertex. Figure 9 shows 7^Q(.(ie„; 0), 
^q=o(ie";0) and T^(ie„;0) defined in eqs. (52) and (54), 
where Q = (tt, • • • , tt) representing the corner of the 
hyper-cubic Brillouin zone. The evaluation involves the 
i-matrices given in eqs. (45) and (47), and nq(ie„; 0) de- 
fined in eq. (59). The imaginary part of nq(ie„; 0) is zero 
due to the particle- hole symmetry at Vm = 0, and the 
real part is shown in Fig. 10. From Fig. 9, the lowest or- 
der functions turn out to be almost identical with each 
other except for low frequencies. Hence, the difference of 
their inverses, Pq = 7^JJ^(ie„; 0)""'^— 7^°(ie„; O)"""^, vanishes 
rapidly with increasing n. This feature of Pq ensures an 
efficient and precise evaluation of the infinite-size matrix 
equation, eq. (56). Only the low- frequency part where Pq 
has finite values is to be provided to obtain the dynam- 
ical susceptibilities of the local moments (see Appendix 
C for detail). This idea is also applicable to evaluations 
of the conduction-electron susceptibility from nq(ie„; 0). 
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temperature ciynamics of the Kondo-type model. A nu- 
merical work in this framework will be presented in an- 
other publication. 
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Appendix A: Derivation of Expression of the 
Generalized t-matrix 

In this Appendix, we derive equations of motion of 
the single-particle and two-particle Green function, and 
demonstrate the expressions of the f-matrix and the gen- 
eralized t-matrix. 

A.l Some useful relations 

We begin with summarizing some relations which are 
necessary in taking derivative involving the time-ordering 
operator. Letting A to D he the fermion or boson, the 
two- and four-operator products are differentiated as 



Fig. 9. Local, uniform (q = 0) and staggered {q = Q) compo- 
nents of r°(ie„; 0) for J = 0.3 and T = 0.02. 



+ 6{n-T2){Tr[A{n),B{n)]±), 



(A-l) 



and 




Fig. 10. Local, uniform and staggered components of llq(i€„;0) 
for J = 0.3 and T = 0.02. 



7. Summary 

We have developed a framework to deal with the dy- 
namics of highly correlated lattice models with localized 
electrons, such as the Kondo and the CS lattice models. 
The t-matrix, which describes the effect of the localized 
moments on the single-particle excitation, has been ex- 
tended to the two-particle responses. The ordinary and 
the generalized t-matrices play essential roles in applying 
the conventional diagrammatical approach to models in 
the atomic limit. Consequently, the spatially dependent 
susceptibilities of the localized moments have been given 
in terms of the generalized t-matrix based on the DMFT. 

Our approach deals with the localized degrees of free- 
dom from the strong correlation limit. The CT-QMC is 
very powerful as an impurity solver which gives finite- 



—{TrAin)B{T2)C{n)D{u)) 

OTl 

T,^^B{r2)C{Ts)D{T,) 

+ S{n - T2){Tr[A{n),B{n)]±C{T3)D{Ti)) 

+ 5{n - Ts){Tr[A{n), C{n)]±D{u)B{T2)) 

+ 6{ri - Ti){Tr[A{Ti), D{ti)]±B{t2)C{t3)), (A-2) 

where [A, B]± = AB ± BA is an anti-commutator and 
a commutator, and arc chosen for fermion and boson, 
respectively. The following relation is satisfied when A 
and B are fermion operators and (p a boson operator: 

A(r,yl(ri)i?(r2)0(r3)) 

a4(Ti) . 

+ S(n-T2){TAA{n),B{Ti)}cf>{T3)) 

+ 5{r, - T3){Tr[A{T,),<P{T,)]B{T2)). (A.3) 

A. 2 Single-particle Green function 

The single-particle Green function is defined by 

Gi2(Ti,r2) = -(T,ci(ri)ct(T2)). (A-4) 

We first differentiate with respect to ti. Using eq. (A-l) 
and substituting eq. (11), we obtain 

d 



Gi2(ti,T2) 

= 5i26{n - T2) - {Trjl{n)cl{T2)). (A-5) 
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The second term is differentiated with respect to T2 to 
yield 



(^-6) {Trh{rM{r2)) 



= {Tr3l{Tl)A{T2)) - S{n - T2)({.n, 4}>- (A-6) 

After the Fourier transform, defined by eq. (10), we ob- 
tain 

Gi2(ie„) = (5i25i(ien) + 5i(ien)ti2(ie„)52(ien), (A-7) 

where the f-matrix is defined by 

ti2{ri,T2) = -(TVii(ri)i|(r2)) +<5(ri -T2)({ii,4}). 

(A-8) 

This equation yields eq. (14) for the impurity model, and 
eq. (16) for the periodic model. 

A., 3 Two-particle Green function 

We proceed to the two-particle Green function in 
eq. (19). As shown in the above derivation, time- 
derivative of c(r) or c^{t) brings about j{T) or j^(t) with 
extra terms. In the following, we in turn differentiate the 
conduction-electron operators until all the operators are 
replaced by j and . 

We first differentiate with respect to ri. Noting the 
relation in eq. (A-2), the equation of motion is given by 

- a) (T.cI(n)c2(r2)4(r3)c4(T4)) 

= (7VjI(Ti)C2(r2)4(T3)C4(T4)) 

+ 5(ri - T2)(5i2G43(t4,T3) - S{ti - T4)(5i4G23(r2, T;,). 

(A-9) 

We introduce an operator which gives the Fourier 
transform as 

jr(T^4(n)C2(T2)4(T3)C4(T4)) 

= ^^ri • • • dr4(r,4(ri)c2(T2)4(r3)c4(T4)) 



(A-IO) 



This transformation replaces the differential operators as 
follows: 

_d_ . d_ . 



(ATI) 



d . . d . 
>ie„/-Fif„, >ie„'. 

Equation (A-9) is then transformed to give 

jr(r^cl(Ti)C2(r2)4(7-3)C4(T4)) 

= gi{\enW{Tr3\{Tl)C2{T2)c\{Ts)Ci{Ti)) 

+ ^mO<5l2G43(ien') — 5nn' 5uG2si}^n + it'm)]- 

(A-12) 



We proceed to make further derivative with respect to 
other time variables T2, ts and T4. We obtain the follow- 
ing equations: 

^(r,iI(Ti)c2(T2)4(r3)c4(T4)) 

= .92(ien + it^m)[-?^(TVjI('^l)j2(r2)4('^3)c4(T4)) 

'523.94(ien)t4i(ie«) 

+ <5mO'53453(ien' + "An) ({.?!, C2})] 

+ ff2(ien + ii^m).9;i(ien' + ii^m),94(ien) 

X {mri - r2){TAjl{ri),C2{ri)}jl{T3)MTi))] 

+ mri - T2)6{n - u){{jl{n),c2{n)}{jl{Ts),c4nm }, 

(A- 13) 

^(r,il(Ti)j2(r2)4(^3)c4(T4)) 

= fl'3(ien' +'^^m){^{Trjl{Tl)j2{T2)jl{T3)C4{T4)) 
+ J^[5{ts - r4)534{Trjl{n)j2{T2))]} 

- fi'3(ien' +it'm)5'4(ien') 

X {^[5(t2 - T3)(T,iI(Ti){j2(r3),4(^3)}i4(T4))] 

+ J'[S{n - T4)S{t2 - r3)(T,{it(Ti),C4(Ti)}{j2(r3),4(^3)})] }, 

(A.14) 

nTrjliri)j2{T2)jl{T3)C4iT4)) 

= 54(ien){^(7ViI (n )j2 (T2)il (^3) j4 (r4)) 

+ mn - r4){TAjKn),C4{r,)}j2{T2)jl{Ts))] 

+ - r4){Trjl{n)j2{T2){jl{rs), C4(r3)})] }. 

(A-15) 

Substituting eqs. (A-13) (A 15) into eq. (AT2) and after 
some manipulations, wc finally obtain an expression of 

J^{Trc\{Ti)c2{T2)cl{Ts)c4{T4)) as follows: 

{Trcl (ti ) C2 (r2 ) 4 (7-3 ) C4 (t4 ) ) 
= '5„io[(5i25i(ieri) • G43(ie„') 

+ .92(ien)^2i(ien).9i(ien) • '534g4(ien')] 
[(51431 (ie«) • G23(ie„ + iz^m) 

+ 54(ien)i4i(ien)5i(ien) • <52352(ien + iz^m)] 
+ fli(ierj)5'2(ien + it'm)5'3(ien' + i^'m)5'4(ien') 

X •^^1234('ri>T-2,T3,r4), (A-16) 

where '7i234(Ti,T2,r3,T4) is defined by 

^1234(n,r2,r3,T4) = (^Vji (Tl)j2(r2)i|(r3)i4(r4)) 

+ S{n-T2)S{T,-T4){Tr{jl{T^),C2{Ti)}{jl{T3),C4{T3)}) 

- 5{t, - T4)5{t2 - r3)(T,0j(ri),C4(ri)}{jt(r3),C2(T3)}) 

+ S{n - T2){Tr{jl{n),C2{ri)}.i{Ts).u{T4)) 

+ S{t3 - r4){Tj,{n)j2{r2){jl{T3), Ciir^)}) 

- 5{n - r4){Tr{jl{n),c4{n)}jl{n)j2{r2)) 

- S{t2 - r3){Trjl{n)j4{T4){jl{Ts), C2(t3)}). (A-IT) 
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In this derivation, we have used {j2, cj} = {02,^3}, which 
is demonstrated from the definition of j2 and j^. Noting 
eq. (A- 7), eq. (A- 16) is rewritten as 

J"(TrCi(Ti)c2(T2)4(T'3)c4(T4)) - (5moG21 (ien)G43 (iCn' ) 
= -Snn'lhigiiien) ' G23(ieri + il^m) 

+ .94(ien)t4i(ie«).9i(ie„) • (523.g2(ie« + i^^m)] 
+ 5i(ie«)52(ien + ij^m).93(ien' + ii'm)54(ien') 

X J^Ti234(Ti,T2,r3,r4), (A-18) 
where Ti234(ti, T2, T3, T4) is given by 

"71234 (ti, r2, T3, T4) 
= ^1234(n,r2,r3,T4) - t21 (t2 , Ti )t43 (t4 , T3) • (A-19) 

As a result, we obtain the expression of 
'7"i234(Ti,T2,r3,r4) in eq. (23). 

Appendix B: Internal Energy 

In this appendix, we derive an expression of the in- 
ternal energy for the impurity and periodic CS models. 
We begin by considering the CS lattice model in eq. (3). 
Although the Hamiltonian include two-body interaction, 
its expectation value can be represented by the single- 
particle Green function with use of the technique intro- 
duced in rcf. 8. Considering an equation of motion for 
Gka,ka{T,T') with t' ^ T + 0, and then taking a sum- 
mation for indices, we obtain 

- lim ■^y^Gka,kaiT,T') 
koL 




dr 



where we have used eq. (11). Taking the Fourier trans- 
form with respect to r, we obtain the internal energy 
E = {HcsL.) as 



Substituting this expression into eq. (B-2), we obtain the 
following equation for a change in the internal energy due 
to the impurity: 

£^imp = (Has) - (-ffc)c 

= ^Ei^" ( l^E^La^n)) ta(ien)e'^"^ (B.4) 

n V ° fca / 



Here, we have set fi = 0. 



of 



High- 



Appendix C: Implementation 
Frequency Limit 

As presented in §4, the q-dcpondcnt dynamical suscep- 
tibility Xq(i^m) in the CS lattice model has been given 
by the high-frequency limit of the generalized t-matrix 
7^(ie„,ie„'; Ii/to). To obtain Tq, we need to solve the 
infinite-size matrix equation, eq. (56). In this appendix, 
we present a efficient way to deal with the infinite-size 
matrices accurately. 

The matrix equation in cq. (56) may be rewritten as 

'^q — %oc + %ocPq%oc + %ocPq%ocPqT\oc H ) (C'l) 

where Pq is a diagonal matrix whose element is defined 
by 

{Pq)i = 7lO,(iQ;ii/„)-i - T,°(iej;ii/„)-i. (C-2) 

In eq. (C-2), only a low-frequency part has finite values, 
since two terms cancel out each other at high frequencies. 
Therefore, we can replace the infinite sum over I with a 
summation up to L 



-t-oo 



(C-3) 



\1\<L 



In Fig. 9, for example, L = 10 is sufficient for reliable 
calculations. With use of this restricted sum, eq. (C-1) is 
rewritten as 

(B-1) {Tq)nn' = {T\oc)nn' + {1\oc)nl{Qq)ll'{T\oc)l'n' , (C-4) 



where the matrix Qq is defined within the restricted 
space as follows: 



^ = T^^ie„Gfc„,fc„(ie„)e'^"*+/i(iVe), (B-2) Qq = Pq+ PqT\ocPq + ■ ■ ■ = - PqT\oc] ^Pq- (C-5) 

Hence, in order to obtain (7q)+oo.-cc. which yields 
Xq(ii^m) in eq. (37), we need to evaluate (Tioc);;' at 



n ka 

where 5 is a positive infinitesimal. In numerical calcu- 
lations, we may subtract l/ie„ from the Green function, 
because it does not contribute to the internal energy. Fur- 
thermore, the quadratic term, c/(ieri)^, can be evaluated 
to be c/2. By expanding eq. (41) in powers of 1/z, we 
obtain c = Jn^ — /i on condition of p{e) = p{—e). Hence, 
excepting this term, the series in eq. (B-2) converges as 
(ie„)~^ at high frequencies. 

Equation (B-2) is also applicable to the impurity CS 
model. In this case, the total energy E = (i?cs) is com- 
posed by the kinetic-energy part being proportional to 
Nq and an impurity contribution. To distinguish the im- 
purity contribution, we express Gfca,fca(ien) in terms of 
the impurity t-matrix ta(ie„) as follows: 

Gfca,fca(ien) = fffc(ien) + fffc(ien)^^^5fc(ien)- (B-3) 

iVo 



low frequencies as well as (7ioc)+ooj, {T\oc)i'-oo, and 
(7ioc)+oo,-oo- The limit e„ ^ oo in TJoc can be taken 
strictly in the CT-QMC. Consequently, the infinite-size 
matrix equation can be solved rigorously to give Xq i}^m ) - 

Appendix D: Vertex Parts for Conduction Elec- 
trons 

We have introduced the generalized t-matrix T from 
the two-particle Green function of conduction electrons 
X^- With use of the integral equation for T, the spatial 
dependence of T in the periodic system has been derived. 
In this appendix, we derive an integral equation for x'^ 
itself. 

The source of the spatial dependence is only ejc 
in the Anderson lattice and the CS lattice models. 
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Hence, we can represent all the spatial dependences ex- 
plicitly by G^(ie„), and therefore by nfeq(ie„; ifm) = 
— G^(ie„)G^^q(ie„ + ivm)- In the following, we omit the 
spin degeneracy for simplicity. The snsccptibility of the 
conduction electrons Xkk'q given in terms of T^k'q in 
eq. (35), which may be rewritten as 

Xfcfc'q(i£nj i^n'; i^^m) = '^rm"5fcfc' Ilfcq (ie„ ; ii^m) 

+ ,9fe(ien),gfc+q(ien + i;^m).gfe'+q(ien' + i^'m)5fe'(ien') 

X ['^fc'q(iern iCn'; it'm) " Snn'Skk'Tkqi^^n'j^'^m)], 

(D-l) 

where T^°q(ie„;iz^TO) = -tk{ien)tk+q{i€n + i'^m)- Using the 
irreducible vertex part /(ie„, ie„/; iz^„i) in cq. (53), the 
terms in the bracket is expressed in the matrix form for 
e„ and e„i as 

%k'q - ^kk'Tkq = '^kql'q'^k'q^ (D-2) 

where /q(ie„, ie„'; ii/m) is a reducible vertex part com- 
posed of I and 7^° : 

/; = J + IT^I + IT^IT^I + ■■■ . (D-3) 

From eqs. (46) and (47), we can derive the relation 
fffc(ien)ife(ie«) = G^(ie„)S'=(ie„), which leads to 



= nfcq(ie„; Wm)T.'' {len)^" {ien + ii 



(D-4) 



Substituting eqs. (D-2) and (D-4) into eq. (D-l), and then 
taking the summation over k and fe', we obtain 

Xq(ien, icn'; ii^m) = nq(ie„; iym)5nn' 

+ nq(ie„; ifTO)ii'q(ie„, ie„'; iz/„)nq(ie„'; ii^m), (D-5) 
where 

-^q(ien, ien';ium) = S'^(ie„)S'^(ie„ + ium) 

X 7^(ie„,ie„';it'„)S^(ie„/)E'=(ie„/ -|-iz/„). (D-6) 

We focus on I'^ defined in cq. (D-3). The spatial depen- 
dence of Ig originates in 7^°, which can be expressed in 
terms of Ilq. With use of the expression of tk in eq. (47), 
defined by eq. (54) is rewritten as 

Tq°(ie„;iz/„) = S=(ie„)2[A(ie„; iz/„) 

+ nq(ie„; ii/„)]S'=(ie„ + iu^f, (D-7) 
where A{ien;ii'm) is defined by 

A{ien; Wm) = 5]'^(ie„)"^S^(ie„ + \Vm)~'^ 
+ S=(ie„)-iG^(ie„ + ii/„) 
+ G^(ie„)S^(ie„ + iz/„)-i. (D-8) 



It is clear from the above equation that A composed only 
of local quantities. Introducing an auxiliary quantity K 

by 

K{ien, ie„'; ii'm) = i;''(ie„)S'=(ie„ -I- ivra) 

X 7(ie„,ie„/;if„)S'=(ie„05^''(ien' +ii'm), (D-9) 
K'q in eq. (D-6) is rewritten as 



K'^ = K + K{A + Ilg)Klj 



(D-10) 

where F'^ is defined by 

= {K-'^ + Ay^ . (D-11) 

Substituting eq. (D-10) into eq. (D-5), we finally obtain 
the equation for Xq as follows: 

x^q = nq + nqF^nq + --- 

= nq + nqF^x',- (D-12) 

Therefore, F'^ signifies the irreducible vertex part for the 
conduction electrons, and is composed only of local quan- 
tities in infinite dimensions. 
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